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1.  INTRODUCTION 


1 . 1  CONSTITUTIVE  BEHAVIOR 

In  the  elastohydrodynamic  (EHD)  regime  of  lubrication,  the  rheology  of  the  liquid 
lubricant  is  key  to  the  generation  of  a  protective  film  and  the  transfer  of  shear  across 
concentrated  contacts  in  many  machine  elements.  It  is  now  understood  that  small  scale 
EHD  plays  a  major  role  in  the  boundary  lubrication  of  rough  surfaces.  Many  high- 
pressure  metal  working  operations  share  the  pressure  and  kinematics  of 
elastohydrodynamics.  A  complete  solution  of  the  EHD  problem,  however,  requires  a 
thorough  knowledge  of  the  lubricant  constitutive  behavior  and  the  attending  properties  as 
functions  of  temperature  and  pressure.  The  Newtonian  assumption  alone  is  often 
inadequate,  and  the  assumptions  upon  which  the  Reynolds  equation  is  based  must  be  re¬ 
examined. 

Early  investigations  of  high-pressure  lubricant  rheology  addressed  the  small  strain 
response.  For  small  strains,  viscous  heating  which  plagues  large  strain  measurements 
may  be  ignored.  These  measurements  were  both  optical  and  acoustical,  e.g.  references 
[1  ]  and  [2]  respectively,  and  provided  linear  viscoelastic  properties.  These  techniques  are 
experimentally  convenient  in  requiring  only  optical  or  electrical  paths  through  the  high- 
pressure  chamber.  The  linear  viscoelastic  response  is,  however,  unable  to  account  for  any 
but  the  small  slip  traction  in  EHD. 

Attention  then  turned  to  large  strain  investigations  which  could  reveal  non-linear 
behavior  and  successfully  predict  EHD  traction  [3].  The  progress  in  this  area  has  recently 
been  significant  in  that  a  mechanism  of  apparent  non-linearity  (at  least  for  low  molecular 
weight  liquids)  has  been  shown  to  be  a  mechanically  induced  shear  localization  or  shear 
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bands  [4],  The  constitutive  behavior  in  the  material  surrounding  the  shear  bands  is 
apparently  linear  viscoelastic. 

Fundamental  to  the  generation  of  a  film  in  concentrated  contacts  is  the 
piezoviscous  property  of  liquid  lubricants  -  whereby  even  simple  low  molecular  weight 
liquids  at  ambient  temperature  attain  very  great  viscosities  under  pressure.  Shear 
rheological  investigations  of  non-Newtonian  response  of  lubricants  [5,6]  have  resorted  to 
low  temperatures  to  achieve  the  same  level  of  viscosity.  It  is  clear  now,  that  high- 
pressure  is  essential  to  an  accurate  simulation.  Liquids  become  considerably  "stronger" 
under  pressure.  Peculiar  to  these  atmospheric  pressure  studies  was  the  observation  of 
fracture.  Eastwood  and  Harrison  [5]  observed  the  liquid  in  shear  and  reported  cracking. 

Early  investigations  of  non-Newtonian  lubricant  response  often  studied  polymer 
solutions.  The  blend  of  mineral  base  oil  and  polymeric  viscosity  index  improver  is 
representative  of  multi-grade  motor  oils.  The  capillary  rheometer  was  useful  for  this 
purpose  as  it  is  very  simply  pressurized  and  although  the  shear  stress  cannot  be  large  for 
this  instrument,  non-Newtonian  flow  occurs  at  relatively  low  stress  for  polymer  solutions. 
The  high-pressure  capillary  viscometer  was  developed  to  a  high  level  by  Jakobsen  and 

Winer  [7]  who  reported  measurements  to  shear  rates  of  10^s'^,  pressures  to  600  MPa  and 
time  of  shear  as  short  as  4  ps.  They  reported  Newtonian  flow  for  liquid  lubricant  base 
stocks  to  a  shear  stress  of  5  MPa. 

1 .2  REYNOLDS  EQUATION 

The  analysis  of  the  pressure  generated  by  lubricant  films  is  almost  exclusively 
performed  with  a  form  of  the  classical  Reynolds  equation.  This  differential  equation 
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derives  from  the  inertialess  form  of  the  Navier-Stokes  equation  combined  with  the 
continuity  equation  with  the  assumption  that  the  flow  channel  is  small  in  one  coordinate 
direction.  It  has  been  generalized  to  incorporate  variable  viscosity  as  well  as  variable 
density  and  has  been  remarkably  accurate  in  predicting  the  film  thickness  in  non- 
conformal  contact  problems.  However,  there  is  a  fundamental  limitation  to  the  Reynolds 
equation  for  problems  in  which  viscosity  varies  with  pressure  and  in  particular,  in  the 
elastohydrodynamic  regime.  The  cross-film  pressure  gradient  cannot  be  neglected  and 
secondary  flows  result. 

Renardy  [8]  recognized  that  the  Navier-Stokes  equations  for  an  incompressible, 
piezoviscous  fluid  may  suffer  from  non-existence  and  nonuniqueness  problems  when  the 
principal  tensorial  strain  rates  are  not  less  than  (2pa)''.  Here,  p  is  viscosity  and 
a  =  d  Inp  /dp  is  the  local  pressure  viscosity  coefficient.  The  Navier-Stokes  equations  can 
undergo  a  change  of  type — a  process  which  has  been  used  to  characterize  shear 
localization  [9]. 

2.  NON-LINEAR  SHEAR  RESPONSE 

The  search  for  the  relevant  constitutive  equations  which  relate  the  lubricant  stress 
to  the  flow  kinematics  in  a  concentrated  contact  has  occupied  the  interest  of  tribologists 
working  in  EHD  for  at  least  thirty  years.  The  goal  has  been  to  construct  experiments  to 
verify  the  rheological  models  and  to  provide  the  necessary  property  relations  so  that  the 
complete  elastohydrodynamic  solution  may  be  obtained.  Various  non-Newtonian  models 
(see  [3]  for  a  review)  have  been  advanced  which  provide  accurate  solutions  for  the 
traction  over  some  operating  range  when  the  required  rheological  properties  are  obtained 
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from  the  same  traction  data.  In  1972,  Dyson  [10]  warned  that  EHD  traction  research  was 
"enclosed  within  a  tight  circle"  of  fitting  parameters  to  observations  without 
consideration  for  measurements  made  outside  of  EHD.  Hopefully,  we  have  broken  from 
the  circle. 

We  may  now  generate  rheological  flow  curves  under  conditions  of  pressure, 
temperature,  and  rate  of  shear  which,  although  still  rather  restricted,  are  sufficiently 
within  the  realm  of  EHD  to  make  accurate  traction  predictions  and  compare  with  traction 
experiments.  We  have  previously  interpreted  these  flow  curves  as  lubricant  constitutive 
behavior.  However,  in  light  of  the  recent  observation  [4]  of  mechanical  shear  bands 
operating  within  the  lubricant  film  concurrent  with  non-linear  shear  response,  this 
interpretation  must  be  accompanied  by  a  caveat:  Rheological  flow  curves  which  are 
generated  in  plane  shear  yield  an  empirical  rate  equation  which  is  useful  in  modeling 
Couette  dominated  lubrication  problems.  A  rigorous  analysis  of  the  EHD  problem  would 
require  a  constitutive  equation  and  a  slip  criterion  such  as  Mohr-Coulomb  [11]. 


2.1  VISCOUS  HEATING  IN  COUETTE  VISCOMETRY 

Presently,  the  most  useful  rheometer  configuration  for  investigation  of  high- 
pressure,  high-stress  response  of  the  liquid  state  is  that  of  rotating  concentric  cylinders. 
When  a  liquid  is  sheared,  the  viscous  work  done  raises  the  temperature  of  the  liquid. 
While  the  study  of  this  phenomenon  is  of  itself  interesting  to  lubrication,  it  is  to  be 
avoided  in  a  rheological  measurement  because  constitutive  behavior  excludes  processes 
which  result  from  temperature  variation.  In  previous  work  [12]  concerning  Couette 
devices  the  authors  have  emphasized  the  importance  of  a  low  Brinkman  No.  (through 
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primarily  a  small  shearing  gap)  to  mitigate  the  effect  of  the  temperature  difference  within 
the  liquid  film  and  fast  instrument  response  to  control  the  temperature  of  the  surfaces  of 
the  solid  boundaries  to  the  film.  An  alternate  experimental  approach  is  to  perform  a 
measurement  in  so  short  a  time  that  the  temperature  profile  in  the  film  has  not  had  time  to 
develop  but  at  a  late  enough  time  that  the  velocity  profile  is  fully  formed.  Winter  [13] 
showed  that  this  latter  approach  requires  that  the  Prandtl  No.  be  large. 

An  earlier  analysis  of  the  combined  effect  of  instrument  response  time  and 
cylinder  heating  [12]  modeled  the  cylinders  as  a  lumped  heat  capacity.  It  was  found  that 
to  minimize  errors  due  to  this  combination  the  measured  stress  history  should  show  no 
thermal  softening  over  a  period  of  time  equal  to  the  instrument  response  time.  Later,  the 
evolution  of  the  radial  temperature  profile  within  the  cylinders  was  considered  [14],  A 
detailed  example  of  a  numerical  thermal  simulation  of  one  of  the  high-pressure  Couette 
viscometers  used  in  this  laboratory  follows. 

Only  conduction  in  the  radial  direction  is  considered.  An  appropriate  form  of  the 
energy  equation  for  the  metal  cylinders  is 

£l  +  l£L,£«&L  £L  (1) 

dr  r  dr  km  dt 

where  r  is  the  radial  coordinate,  t  is  time  and  p,  C,  and  k  are  the  density,  specific  heat 
capacity  and  thermal  conductivity,  respectively.  The  subscript,  m,  refers  to  the  metal 
from  which  the  cylinders  are  fabricated.  Within  a  pressure  vessel  the  outer  surface  of  the 
outer  cylinder  (or  cup)  is  in  contact  with  low  conductivity  liquid.  Therefore  at  the 
outermost  surface  the  boundary  condition  is  set  at  dT/dr  =  0,  namely  it  is  an  adiabatic 
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surface.  Symmetry  at  that  axis  of  rotation  dictates  a  similar  boundary  condition,  cT/dr  = 
0  at  r  =  0. 

The  shearing  gap  is  very  thin  compared  with  the  working  radius  so  that  curvature 
of  the  liquid  film  can  be  ignored.  A  preliminary  numerical  solution  for  the  liquid  film 
with  isothermal  boundaries  and  dissipation  showed  that  for  a  1  pm  film,  steady  state 
temperature  distribution  was  achieved  after  10  ps  of  shearing.  Therefore,  the  storage 
term  in  the  energy  equation  for  the  liquid  film  is  dropped.  Including  the  dissipation  term, 
tv  .  we  have  for  the  liquid 

+  TY  *  0  (2) 

or 

The  geometry  considered  is  shown  in  Figure  1 .  The  liquid  being  sheared  is  in  the 
circular  gap  between  the  outer  and  inner  cylinders.  Only  the  outer  cylinder  rotates  with  a 
velocity  which  yields  a  specified  average  of  the  rate  of  shear,  y .  The  temperature  and  the 
radial  heat  flux  are  made  continuous  at  the  two  metal/liquid  interfaces.  Equations  (1 )  and 
(2)  in  their  respective  regions  are  solved  numerically  with  the  above  boundary  conditions 
and  with  initial  uniform  temperature  (arbitrarily  set  to  zero  for  Figure  1)  to  obtain 
temperature  distribution  and  shear  stress  versus  time. 

The  result  of  a  sample  computation  is  shown  in  Figure  1  for  two  velocity 
histories.  The  velocity  of  the  outer  cylinder  was  increased  linearly  with  time  until  an 

apparent  or  average  rate  of  shear  of  lO^s’l  was  reached  at  to  equal  to  10  and  50  ms.  The 

liquid  rheology  was  Newtonian  with  p  =  10^  Pa-s  and  p  =  0.1  °C*1.  The  gap  was  1.0  pm. 
The  liquid  and  metal  thermal  conductivities  were  k  =  0.15  and  km  =  225  W/m°C  and  for 

the  metal  pm  Cm  -  3  x  106  J/m^°C.  The  temperature  profiles  within  the  metal  cylinders 
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are  shown  in  Figure  1  for  t  =  t@.  The  liquid  film  temperature  reached  a  maximum  very 

near  the  surface  of  the  inner  cylinder.  The  reduction  in  shear  stress  from  isothermal  was 
less  than  2%  for  t0  =  10  ms  and  less  than  4%  for  t0  =  50  ms.  The  rotational  velocity  in 

this  example  is  only  about  1  revolution  per  second  and  can  easily  be  achieved  in  10  ms  by 
a  stepper  motor.  When  a  very  large  driving  torque  is  required  a  50  ms  linear  velocity 
ramp  may  be  necessary. 

In  early  high-pressure  Couette  devices  the  critical  response  time  was  often  that  of 
the  torque  transducer  which  provides  the  stress  measurement.  This  transducer  is 
immersed  in  a  high-pressure  liquid  and  for  high  viscosity  the  transducer  responds  as  a 
first  order  instrument  w-ith  a  time  constant,  tc,  which  is  proportional  to  liquid  viscosity. 
Very  fast  measurements  may  be  obtained  by  accelerating  both  cylinders  and  then 
arresting  the  rotation  of  one  cylinder  with  the  torque  transducer  which  can  be  immersed 
in  a  low  viscosity  liquid.  Now',  the  velocity  response  is  intimately  connected  with  the 
transducer  response.  The  undamped  natural  frequency  for  the  transducer  is  2.5  kHz. 
Theoretically,  a  5  percent  settling  time  [15]  of  0.2  ms  can  be  achieved  for  transducer 
response  by  optimizing  the  damping  liquid  viscosity.  This  is  difficult  to  achieve  in 
practice  as  the  viscosity  of  the  damping  liquid  changes  with  the  test  temperature  and 
pressure.  In  any  event,  we  have  not  observed  an  improvement  in  measurements  typical  of 
those  reported  here  when  done  faster  than  50  ms.  The  temperature  increases  shown  in 
Figure  1  present  no  problem  in  interpretation  of  results. 

Also,  by  not  considering  axial  conduction  in  the  cylinders  the  analysis 
overestimates  the  metal  temperature  rise.  Comparing  measured  stress  histories  for  a 
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Newtonian  liquid  with  the  computed  stress  history  at  t  =  Is  we  found  the  actual  loss  of 
shear  stress  compared  to  isothermal  shear  was  only  64%  of  the  predicted  loss. 

The  temperature  distribution  of  Figure  1  was  incorporated  into  a  thermal-elastic 
numerical  analysis  of  each  cylinder  using  both  plane-strain  and  plane-stress  idealizations. 
Thermal  expansion  results  in  no  subtantial  change  in  the  operating  film  thickness  at  the 
time  of  the  measurement. 

2.2  SOME  RESULTS  FOR  PRESSURES  TO  300  MPA 

A  high-pressure,  high  shear  stress  Couette  viscometer  was  developed  [14]  in  this 
laboratory  to  operate  at  pressures  to  300  MPa.  The  steady  shear  response  of  various 
liquid  lubricants  [14]  and  a  grease  [16]  have  been  investigated.  For  the  grease,  an 
empirical  stress  equation  was  found  of  the  form 

Tq  =  Tj  +  (gpy)q  ,T  >  (3) 

where  g  is  a  dimensionless  viscosity  enhancement  factor  approximately  equal  to  2,  q  is  a 
dimensionless  exponent  and  p  is  the  base  oil  viscosity.  For  q  =  0.5,  equation  (3)  becomes 
Casson's  equation  for  oil/particulate  mixtures. 

Flow  curves  are  presented  for  a  high-traction  fluid,  Santotrac  50  and  a  mineral  oil, 
LVI  260  in  Figures  2  and  3,  respectively.  Although  the  pressures  and  temperatures  are 
comparable  the  shear  response  is  quite  different.  The  cycloaliphatic  traction  fluid 
remains  Newtonian  to  a  shear  stress  of  at  least  10  MPa  whereas  the  mineral  oil  responds 
non-linearly  above  3  MPa.  Note  also  that  for  the  mineral  oil,  the  stress  becomes 
essentially  independent  of  rate  at  about  x  =  10  MPa  and  that  the  temperature  dependence 
diminishes  as  the  rate  dependence  diminishes.  Because  of  their  extended  Newtonian 
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response  we  have  begun  to  use  the  cycloaliphatic  traction  fluids  as  viscosity  standards  for 
shear  stress  greater  than  1  MPa  in  calibrating  cylinder  gaps. 

2.3  A  NEW  HIGH-STRESS  COUETTE  VISCOMETER  FOR  600  MPA  PRESSURE 

An  increase  in  the  pressure  capability  of  previous  rotating  concentric  cylinder 
viscometers  is  necessary  to  investigate  the  shear  response  of  typical  lubricants  well  above 
the  glass  transition  temperature.  Translating  concentric  cylinder  rheometers  are  capable 
of  pressures  above  1  GPa;  however  the  rate  of  shear  is  so  low  that  the  viscosity  must  be 
made  very  large  in  order  to  achieve  an  interesting  magnitude  of  shear  stress.  Hence 
experiments  are  often  conducted  near  the  glass  transition  of  the  lubricant. 

We  recently  reported  [14]  the  development  of  a  Couette  viscometer  for  pressure  to  300 
MPa  and  the  results  of  the  previous  section  were  obtained  with  that  device.  The  pressure 
capability  of  the  Couette  technique  is  doubled  with  the  device  shown  in  section  view  in 
Figure  4.  The  outer  cylinder  (cup)  is  driven  by  an  external  stepper  motor  by  means  of  a 
drive  shaft.  The  thrust  bearing  which  prevents  expulsion  of  the  shaft  from  the  vessel  has 
been  moved  to  outside  of  the  vessel  so  that  the  bearing  does  not  run  in  pressurized  liquid. 
The  high-pressure  seal  is  now  a  simple  spring-loaded  packing.  There  is  no  need  for  high- 
pressure  electrical  feed-throughs  since  the  torque  measurement  is  transmitted  out  of  the 
vessel  optically.  The  inner  cylinder  (bob)  is  restrained  from  rotation  by  a  torsion  bar  with 
a  mirrored  surface  within  a  glass  tube.  See  Figure  4. 

The  principle  of  the  torque  measurement  is  depicted  in  Figure  5.  The  twist  of  the 
torsion  bar  results  in  a  deflection  of  a  laser  beam  which  is  detected  by  an  optical  position 
sensor.  A  glass  tube  (Figure  4)  provides  a  circular  interface  between  the  liquid  which 
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surrounds  the  mirror  and  the  pressurizing  fluid  so  that  changes  in  refractive  index  of  the 
sample  has  no  effect  on  the  measurement. 

2.4  RESULTS  FOR  PRESSURES  TO  600  MPA 

A  flow  chart  is  presented  for  the  polyalphaolefin,  SHF  1001,  in  Figure  6.  The 
response  is  apparently  Newtonian  to  about  2  MPa  shear  stress,  although  the  effective 
viscosity  below  2  MPa  is  slightly  reduced  compared  to  falling  body  measurements. 

The  mineral  oil,  HVI  650,  has  been  the  subject  of  numerous  traction  experiments. 
The  shear  response  of  HVI  650  is  depicted  in  Figures  7  a  and  b.  The  curves  drawn 
through  the  data  points  for  23°C  and  44°C  represent  the  empirical  Carreau-Yasuda 
equation 

t  =  py 

Here,  a  is  a  dimensionless  parameter  which  controls  the  breadth  of  the  transition  from 
Newtonian  to  a  non-Newtonian  flow  regime  with  rate  sensitivity  coefficient  of  m.  For  m 
=  0,  tL  is  a  limiting  shear  stress.  In  Figures  7  a  and  b,  we  used  a  =  1.5  at  23°C  and  a  =  2 
at  44°C.  For  the  limiting  stress,  the  form  proposed  by  Bezot,  et  al.  [17]  was  adopted 

tl  =  c0  +  c,p  +  c2p2  (5) 

Using  the  results  of  Figure  7b  and  Ref  [18]  and  [19],  we  obtained  ci  =  0.034  and  C2  =  2.1 
x  10‘5  MPa  with  c0  arbitrarily  set  to  zero.  Note  that  for  HVI  650  at  44°C  (Figure  7a)  the 
curves  are  also  a  good  approximation  of  the  Sinh  Law, 

t  =  t0  sinh‘](py /t0)  (6) 


1  + 


f  ■  y 

py 

Ul; 


\  jii-i  )■  a 


(4) 
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where  t0  is  a  stress  which  represents  the  limit  of  Newtonian  response  and  is  3.3,  5.2  and 
6.5  MPa  for  pressures  of  448, 5 1 7  and  586  MPa  respectively.  That  is,  x0  *  xL  /  4  at  44°C. 

The  shear  response  of  the  traction  fluid,  Santotrac  40,  is  shown  in  Figure  8.  For 
the  pressure  of  414  MPa,  the  solid  curve  is  equation  (4)  and  the  two  broken  curves  are 
equation  (6)  for  x0  of  6  and  25  MPa  respectively.  When  elastohydrodynamic  traction 
curves  are  interpreted  in  terms  of  equation  (6)  together  with  exponential  pressure- 
viscosity  the  traction  gradient  can  be  set  equal  to  x0  [18] 


5r 

dtny 


=  T„ 


(7) 


Values  of  t0  obtained  from  traction  tests  using  (7)  are  typically  less  than  6  MPa 
[18,  20]  for  cycloaliphatic  traction  fluids  under  conditions  of  Figure  8.  Clearly,  from 
Figure  8.  the  Newtonian  limit  exceeds  6  MPa  and  is  close  to  25  MPa.  Therefore,  any 
interpretation  of  the  traction  gradient  as  a  Newtonian  limit  through  a  Sinh  Law  model 
with  t„  independent  of  pressure  must  be  suspect.  However,  if  the  pressure  of  a 
rheological  measurement  is  chosen  carefully  it  should  be  possible  to  obtain  agreement 
between  the  Newtonian  limit  for  that  particular  pressure  and  the  traction  gradient  since 
the  Newtonian  limit  is  pressure  dependent  (eg.,  HVI 650  at  448  MPa  and  44°C). 

We  have  investigated  the  grease  from  Ref  [16]  and  its  base  oil  in  the  current 
rheometer  to  explore  departures  from  Newtonian  behavior  at  high  shear  stress.  In  Figure 
9,  for  room  temperature,  the  mineral  base  oil,  600P,  is  non-Newtonian  above  t  =  2  MPa 
at  3 1 0  MPa  pressure.  When  the  pressure  is  increased  to  517  MPa  this  oil  is  Newtonian  to 
t  =  3  MPa.  The  CA7000,  which  is  a  soap-thickened  grease  of  600P,  is  also  shown  in 
Figure  9.  It  would  appear  from  these  limited  data  that  the  shear  stress  for  the  grease  (at 
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high  stress)  may  be  obtained  by  multiplying  the  base  oil  result  by  the  previously  [16] 
obtained  viscosity  enhancement  factor,  g,  which  is  about  2.  Clearly,  equation  (3)  is  not 
applicable  here  since  it  becomes  approximately  Newtonian  at  high  shear  stress.  Other 
techniques  have  been  developed  to  probe  the  response  of  liquids  to  high  pressure  and 
stress. 

Other  laboratories  have  observed  rate  independence  at  elevated  pressure.  For 
example,  the  High-Pressure  Impact  Viscometer  of  Wong,  et.  al.  [21]  entraps  a  quantity  of 
liquid  between  a  ball  and  a  plate.  Interferometry  is  used  to  determine  the  local  flow  rate 
of  liquid  leaking  from  the  entrapment  and  local  surface  distortion  which  yield  local 
pressure  from  elasticity  theory.  A  Rabinowitsch  Correction  for  slit  flow  gives  the  true 
shear  rate  at  the  surface.  The  authors  concluded  that  the  effective  viscosity  was  found  in 
earlier  works  to  be  a  function  of  time  because  the  shear  rate  varied  with  time  while  the 
shear  stress  remained  at  the  limiting  value.  Limiting  stress  type  behavior  was  observed 
over  4  decades  of  shear  rate  with  m  *  0.01  for  LVI 260. 

3.  EHD  TRACTION  CALCULATION 

We  should  expect  our  property  measurements  to  yield  reasonable  predictions  of 
concentrated  contact  traction.  Evans  [22]  generated  isothermal,  line-contact  traction 
curv  es  for  three  of  the  most  widely  investigated  liquid  lubricants:  5P4E,  Santotrac  50  and 
HVI650.  These  data  were  described  as  isothermal  since  the  disc  temperature  was 
adjusted  to  provide  a  constant  estimated  average  film  temperature.  These  traction  data 
are  presented  for  HVI650  in  Figure  10  for  inlet  temperatures  of  40  and  60°C  and  average 
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pressures,  p ,  of  0.47  and  0.63  GPa.  The  sliding  velocity  is  AU  and  inlet  temperature  rise 
is  3°C. 

For  the  traction  calculations  which  are  the  curves  in  Figure  10,  we  assumed  the 
Hertzian  pressure  distribution  and  integrated  the  Carreau-Yasuda  equation  (4)  across  the 
contact  area  to  obtain  the  average  shear  stress,  7 .  The  dimensionless  parameter,  a,  was 
specified  by  a  =  63/(T+19°C)  which  yields  the  values  obtained  in  the  previous  section 
without  approaching  the  meaningless  condition  of  a  =  0  at  ordinary  temperatures. 
Viscosity  was  obtained  from  a  Free  Volume  Model  [23]  and  limiting  stress  was  obtained 
from  equation  (5).  Results  were  insensitive  to  the  selection  of  the  rate  sensitivity,  m, 
from  0  to  0.03.  Although  the  predicted  curve  at  40°C  rolls  over  more  quickly  than  the 
measured  traction,  the  general  agreement  is  good.  Considering  the  great  differences 
between  the  two  techniques  and  the  assumptions  involved  (eg.,  dry  Hertzian  pressure 
distribution)  it  might  be  unreasonable  to  expect  to  do  better. 
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4.  MECHANICAL  SHEAR  BANDS  IN  NON-VISCOMETRIC  FLOW 


It  has  long  been  suspected  that  the  rate  independent  behavior  observed  and  discussed 
above  was  related  to  slip;  either  at  the  boundaries,  or  internal  in  the  material.  Visualization  of 
shear  bands  has  to  date  been  accomplished  between  parallel  plates  [4].  The  interpretation  of  the 
flow  field  and  stress  state  leading  up  to  localization  is  the  least  ambiguous  for  plane  Couette 
shear.  This  is  not  however  the  general  case  for  lubrication  flows.  To  investigate  the  generation 
of  shear  bands  in  wedge  flow  we  fabricated  a  new  stationary  shaft  for  the  High-Pressure  Flow 
Visualization  Cell  [4].  The  stationary  surface  (the  lower  surface  of  Figure  11)  was  ground  at  an 
angle  of  5°  with  the  direction  of  motion  of  the  moving  (upper  of  Figure  1 1)  surface.  Observ  ation 
of  birefringence  during  shear  through  crossed  polarizers  clearly  showed  a  shear  stress  gradient 
across  the  film  as  expected.  Shear  bands  appeared  as  shown  in  Figure  1 1  for  5P4E.  These  bands 
are  similar  to  those  observed  between  parallel  plates.  However,  the  Mohr-Coulomb  analysis  can 
not  easily  be  applied  because  of  the  non-uniform  stress  field  in  the  film.  In  solid  mechanics 
isothermal  shear  bands  are  often  associated  with  a  change  in  type  of  the  governing  partial 
differentia]  equations.  Singularities  in  the  governing  equations  accompany  this  change  of  type. 
As  shown  below,  singularities  do  occur  in  the  Navier-Stokes  equation  with  realistic  lubricant 
properties. 


5.  A  FUNDAMENTAL  LIMITATION  OF  THE  REYNOLDS  EQUATION 

FOR  PIEZO  VISCOUS  LIQUIDS 

We  have  just  seen  that  shear  bands  will  develop  during  non-Viscometric  shear  of  a 
lubricant  at  high  pressure  and  high-shear  stress.  These  visible  bands  effectively  represent  slip 
planes  and  are  not  predicted  by  current  numerical  analyses  of  lubricant  behavior. 

The  analysis  of  the  pressure  generated  by  lubricant  films  is  almost  exclusively  performed 
with  a  form  of  the  classical  Reynolds  equation.  This  differential  equation  derives  from  the 
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inertialess  form  of  the  Navier- Stokes  equation  combined  with  the  continuity  equation  with  the 
assumption  that  the  flow  channel  is  small  in  one  coordinate  direction.  It  has  been  generalized  to 
incorporate  variable  viscosity  as  well  as  variable  density  and  has  been  remarkably  accurate  in 
predicting  the  film  thickness  in  non-conformal  contact  problems.  However,  there  is  a 
fundamental  limitation  to  the  Reynolds  equation  for  problems  in  which  viscosity  varies  with 
pressure  and  in  particular  the  regime  of  elastohydrodynamics.  The  cross-film  pressure  gradient 
cannot  be  neglected,  and  secondary  flows  result  even  for  flow  between  parallel  plates. 

Renardy  recognized  [8]  that  the  Navier-Stokes  equations  for  an  incompressible, 
piezoviscous  fluid  may  suffer  from  non-existence  and  nonuniqueness  problems  when  the 
principal  tensorial  strain  rates  are  not  less  than  (2 pa)'1 .  Here,  p  is  viscosity  and  a=d  In  p/dp  is 
the  local  pressure  viscosity  coefficient  where  p  is  the  pressure.  The  Navier-Stokes  equations  can 
undergo  a  change  of  type  -  a  process  which  has  been  used  to  characterize  shear  localization  [9]. 

Recently,  Bair  and  Khonsari  [24]  reported  the  occurrence  of  singular  pressure  gradients  in 
two-dimensional,  inhomogeneous  flows  of  incompressible  piezoviscous  liquids.  In  the  case  of 
flow  between  parallel  plates,  these  singularities  were  shown  to  take  place  when  the  shear  stress 
t  approaches  1/a.  This  criterion  is  equivalent  to  that  found  by  Renardy.  Note  that  the  one¬ 
dimensional  Reynolds  equation  does  not  yield  this  singularity.  In  fact  it  predicts  a  trivial  solution 
for  the  pressure  if  pressure  boundary  conditions  are  fixed  at  a  constant  value.  We  will  show  that 
the  pressure  dependence  of  viscosity  leads  to  very  large  cross  film  pressure  gradients  and 
secondary  flows. 


5.1  PIEZOVISCOUS  REWOLDS  EQUATION 

Consider  the  derivation  of  a  generalized  Reynolds  equation  by  Dowson  [25].  The 
coordinate  system  is  defined  so  that  z  is  in  the  direction  of  the  film  thickness  h.  The  x  and  y 
dimensions  of  the  bearing  are  of  the  order  of  L  which  is  several  orders  of  magnitude  greater  than 
h.  Following  Dowson,  the  first  stage  of  the  order  of  magnitude  analysis  of  the  Navier-Stokes 
equations  for  a  Newtonian  liquid  yields: 


dp__d_  du 
fo  =  5z  dz} 


(8) 


and  a  similar  equation  for  dp/dy. 
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For  the  cross-film  direction,  the  appropriate  equation  is: 


dp  d  du 

dz  ~“'+dx  ^dz 


(9) 


where  we  will  assume  that  the  terms  represented  by  ...  are  of  lesser  magnitude  than  the  term 
retained.  In  arriving  at  these  equations  it  is  apparently  assumed  that  in  general  a  function  such  as 
d(p/d£,  is  of  the  order  oftp/E,  Then  if  z  is  of  the  order  ofh0,  comparing  only  the  right  hand  sides  of 
(8)  and  (9)  gives  the  result  that  (dp/dz)/(dp/dx)  is  of  the  order  of  h0/L.  This  result  is  usually  used 
to  justify  the  omission  of  dp/dz.,  but  clearly  the  left-hand  side  gives  the  opposite  result. 

To  illustrate  the  significance  of  the  cross  flow  pressure  gradient,  we  carry  out  the 
differentiation  of  the  products  in  the  right-hand  sides  of  (8)  and  (9)  as  shown  below: 


dp  dp  du 

—  =  ap— —  + 
dx  dz  dz 


dp  dp  du  d2u 

—  =  — h  au— — +  p— — 
dz  dx  dz  dxdz 


(8a) 

(9a) 


Begin  with  an  inspection  of  (9a).  We  neglect  the  term  involving  cross  derivatives  since  du/dx  is 
expected  to  be  small;  however,  our  argument  is  not  compromised  if  this  term  is  included. 
Applying  a  similar  order-of-magnitude  rule  as  before  shows  that  the  ratio  (dp  /  dz)  /  (dp  /  dx)  is  of 

the  order  of  apU  /  h0.  A  conservative  value  of  a  is  10‘Pa  '.  Taking  the  velocity  in  the  x- 
direction,  U,  to  be  1  m/s  and  hc  to  be  10 7  m,  we  find  that  with  a  viscosity  evaluated  at  low 
pressures,  say  p=10':  Pa.s,  the  order  of  magnitude  of  (dp  /  dz)  /  (dp  /  dx)  is  10'\  This  is  typical  of 
hydrodynamic  lubrication  or  EHL  inlet  zones  and  this  order  of  magnitude  analysis  can  be  used  to 
neglect  the  pressure  gradient  in  the  cross  film  direction.  For  pressures  relevant  to  Hertzian  zone, 
p=104  Pa.s  is  not  unusual  and  the  order  of  magnitude  of  (dp  /  dz)  /  (dp  /  dx)  becomes  10'. 
Clearly,  therefore,  at  high  pressures  dp  /  dz  is  not  insignificant  compared  to  dp  /  dx . 

In  deriving  the  classical  Reynolds  equation  based  on  the  reduced  form  of  the  Navier- 
Stokes  equation,  it  is  typical  to  define  a  dimensionless  viscosity  as  p  =  p/p0,  where  p  is 
assumed  to  be  of  the  order  of  one.  If  p0  is  taken  to  be  ambient  viscosity,  then  p  can  be  very 
large  at  high  pressures.  Even  when  the  cross  film  gradient  is  significant,  it  must  have  an  affect  on 
the  pressure  profile  calculation  for  us  to  question  the  use  of  Reynolds  equation  for 
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elastohydrodynamic  applications.  Notice  that  in  equation  (8a)  the  term  dp  /  dz  appears  explicitly 
and  that  it  is  coupled  to  dp  /  dx .  This  coupling  can  result  in  a  singularity  in  the  pressure  gradient 
which  the  Reynolds  equation  is  incapable  of  predicting. 


5.2.  LUBRICANT  PROPERTIES 

Before  developing  some  properties  of  the  isothermal  piezoviscous  Navier-Stokes 
equations  for  the  Hertzian  zone,  it  is  instructive  to  examine  the  viscosity  and  compressibility  of 
liquid  lubricants  at  very  high  pressure.  Free  volume  models  are  useful  for  lubrication  analyses 
because  they  provide  a  link  between  volume  and  viscosity.  In  the  following  example,  we  make 
use  of  the  isothermal  free  volume  model  of  Cook,  et  al.  [27]. 

The  Dolittle  equation  relates  viscosity  to  volume,  V,  and  occupied  volume,  V^.. 


V- V 


BV„ 


V  -V 

*  0  occ 


+  In  Pc 


(10) 


The  Tait  equation  provides  an  expression  of  the  volume  variation  with  pressure: 


V  1  1  . 

—  =  1 - In 

v  k:  +  i 


'i+Pt!±£! 

K„ 


(ID 


V  J 

At  zero  pressure  the  volume  is  V0,  the  bulk  modulus  is  K0  and  the  pressure  rate  of  change  of  the 
bulk  modulus  is  K.;, .  The  secant  compressibility  is  then 


1 


(12) 


v0dP  k0+(i+k;)P 

It  can  be  shown  empirically  that  Vocc  is  nearly  independent  of  pressure.  Differentiating  (10) 
yields: 


V 

R— 2« 

v 

a  = - £ - 7<J>$ 

f  V  V  ^ 

v  OCC 

vV0 "  V0  , 


(13) 


Parameters  typical  of  a  mineral  oil  at  a  temperature  of  75°C,  are  B=3,  K0=1.5  GPa  K’ =10  and 

Voc/Vo=0.70. 

A  plot  of  <I>S  and  a  as  a  function  of  pressure  is  shown  in  Figure  12.  First  note  the 
remarkable  reduction  of  compressibility  for  increasing  pressures.  Although  the  existence  of  high 


17 


pressure  is  often  a  justification  for  a  compressible  solution,  if  we  restrict  our  attention  to 
pressures  greater  than  1  GPa,  an  incompressible  analysis  is  acceptable.  Note  also  that  following  a 
slight  decrease  at  low  pressures,  a  increases  with  pressure  for  pressures  greater  than  about  0.5 
GPa.  This  is  a  result  of  the  competition  between  the  rapidly  decreasing  compressibility  at  low 
pressures  which  decreases  the  sensitivity  of  viscosity  to  pressure  and  the  reduction  of  free 
volume  at  high  pressures  which  increases  the  sensitivity  of  viscosity  to  pressure.  (Liquids  for 
which  K0  is  particularly  high  will  not  show  the  initial  decrease  of  a  with  p.)  The  pressure- 
viscosity  coefficient  can  be  very  large  at  high  pressures.  This  is  in  contradiction  with  empirical 
rules  often  used  in  EHD  numerical  simulations  for  which  a  is  allowed  to  monotomically 
decrease  with  pressure. 

The  analyses  which  follow  pertain  to  applications  where  the  pressure  is  high  enough  so 
that  compressibility  of  the  liquid  can  be  neglected  and  that  the  shear  stress  is  comparable  to  the 
reciprocal  of  a.  To  make  the  piezoviscous  Navier- Stokes  equations  tractable,  we  will  treat  a  as  a 
constant,  recognizing  that  this  assumption  is  an  idealization. 
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5.3  Navier-Stokes  Equations  for  Two-Dimensional,  Intertialess  Flow  of  an  Incompressible, 
Piezoviscous  Liquid 

Consider  the  creeping  flow  of  a  liquid  for  which  the  flow  channel  is  very  large  in  one  direction, 
y.  The  component  of  the  velocity  in  the  y-direction,  v,  and  all  derivatives  with  respect  to  y  must 
be  zero  to  provide  symmetry.  The  appropriate  Navier-Stokes  equations  are: 


dp  d 
dx  dx 

dp  _  d 
dz  dz 


du  (  2  Y du  dw') 


2^+ 


du 


3  )\dx  dz )]  dz 


2fS  +  |S--M 


^du  dw 
\dx  dz , 


( dw  du') 


(  dw  +  du 


dx|_  v^x  dzJJ 


(14) 

(15) 


where  q  is  the  bulk  viscosity  which  is  not  zero  for  liquids. 

The  continuity  equation  for  an  incompressible  liquid  reads: 

du  dw 


■  *4 - =  0. 

dx  dz 


(16) 


Using  equation  (9),  the  Navier-Stokes  equations  (7)  and  (8)  reduce  to  the  following  form: 


dx  &xyr  x/  dz 1 

where  subscripts  on  velocities  represent  partial  differentiation. 
Note  that 

dp  dp  ,  dp  dp 
dx  ^  dx  dz  dz 

Substitution  and  differentiation  gives 

^  =  2pcx  ux  +  pa(uz  +  +  pV2u 


(17) 

(18) 


dp_„ _  pP 

dz 

Components  of  the  strain-rate  tensor  are: 


^  =  2pct  wz  pa(u2  +  wx)^  +  pV2w 


(19) 

(20) 
(21) 


du; 

dXj  dXj 


(22) 


with  principal  values  of  0,  d„  and  -d„  where  for  this  problem 
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(23) 


-d  d 

XX  zz 


Equations  (20)  and  (21)  can  be  combined  to  give  the  following  result  which  uses  the  above 


notation. 

dp  p[2  p  a  d^V^  +  (1  -  2  p  a  da  )V2u] 
dx  l-(2pad,)2 

dp_  p[2padxzV2u  +  (l-2padx;i)V2w] 
dz  l-(2pad,)2 


(24) 

(25) 


Clearly  for  a  =  0,  —  =  pV2u  and  —  =  pV2w  as  required.  Note  that  there  may  be  singular 
dx  dz 


pressure  gradients  for 


2pa 


(26) 


An  equivalent  expression  for  (26)  is: 

1 


where  t,  is  the  maximum  principal  shear  stress. 


(27) 


5.4  Example 

Consider  the  two-dimensional  flow  of  an  incompressible  piezoviscous  liquid  between 
two  parallel  plates  with  in-plane  realtive  motion.  The  cross-film  direction  is  the  z-direction  with 
thickness  h,  and  plate  motion  is  in  the  x-direction.  For  the  isoviscous  case  the  velocity  profile  is 
parabolic  and  the  x  component  of  velocity  is  a  function  of  z  only.  For  no  wall  slip  and  z=0  at  one 
plate  and  for  an  isoviscous  liquid,  the  components  of  the  velocity  should  be  of  the  forms  given 
below: 

u(z)  =  az  +  bz2 

w  =  0.  (28) 

It  has  been  customary  to  assume  the  velocity  field  given  by  (28)  for  piezoviscous  fluids.  We  will 
show  later  that  in  the  case  of  flow  between  parallel  plates  a  secondary  flow  must  exist  when  the 
viscosity  is  a  function  of  pressure.  For  illustrative  purposes,  however,  we  will  assume  that  w  s  0 
is  a  good  approximation. 
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Equation  (24)  and  (25)  can  be  written  in  terms  of  velocities  as: 


=  D  K1 "  2  a  ^  W*)y2u  +  a  ^  (w* +  ^2w 

|  (24) 

*  ‘^=p’[(1-2aMu*)v2w  +  aM(wx  +  uz)^2u] 

(25) 

with 

D  =  l  +  (pa)2[4uxw2  -(w,  +uz)2] 

(29) 

Substituting  w  =  0,  we  obtain 

dp  M«n 

dx  l-(apu2)2 

(30) 

ap  ap2  U2 ua  ap 

—  =  — - — — sV  =  apu2-^- 
8z  l-(apu2)*  dx 

(31) 

Introduce  the  dimensionless  quantities 

P  =  (p-pR)a:  x  =  (x-xr)pR  u22a;  p  =  — =  ep;  T 

=  u2apR,  (32) 

Mr 


where  subscript  R  denotes  a  reference  value  at  x  =  0. 

We  will  solve  equation  (30)  along  a  constant  z  plane.  Using  the  dimensionless  quantities  defined 
in  (32).  equation  (30)  reads: 

4  = - ^  (30a) 

ax  i-tVp 

which  is  separable.  Integrating  (30a)  and  making  use  of  the  boundary  condition  p(x  =  0)  =  0, 
leads  to 

x  =  l  +  T2-e-p-TV  (33) 

which  can  be  solved  for  p  as 

2T 

where  S  =  x  - 1  +  T2 . 

Figure  13  depicts  the  variation  of  p  as  a  function  of  x  for  T  »  0  and  T  =  0.5.  Notice  that 
for  x  >  (1  -  T)2 ,  there  is  no  real  solution  for  the  pressure.  Notice  also  that  there  are  two  branches 
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of  p(x),  although  the  upper  branch  may  be  discarded  as  it  does  not  satisfy  the  boundary 
condition. 

Turning  our  attention  now  to  the  Reynolds  equation  for  the  piezoviscous  flow  of  a  liquid 
between  parallel  plates  we  can  write 


d  fj_dp 
p  dx 


dx 


Integration  of  (35)  results  in  the  following. 


=  0 


~  =  Cp  =  pu2 
dx 


(35) 


(36) 


where  the  constant  of  integration,  C,  is  un.  Therefore  in  dimensionless  form,  equation  (36) 
reads: 


*E  =  e'> 


dx 


(36a) 


The  above  solution  is  the  same  as  that  predicted  by  equation  (30a)  with  T=0,  which  is  shown  in 
Figure  13.  Comparison  of  the  right-hand-side  of  (30)  with  that  of  (36)  reveals  that  the  Reynolds 
equation  adequately  captures  the  mechanics  of  a  piezoviscous  liquid  only  when  apu2  «  1  (in 
other  words  a  x  «  1 )  assuming  w=0. 


5.5.  SECONDARY  FLOWS 

In  the  previous  example,  it  was  tacitly  assumed  that  for  flow  between  parallel  plates,  the 
component  of  the  velocity  normal  to  the  boundary,  w,  is  zero.  This  is  true  for  isoviscous  flows 
and  can  be  easily  verified.  For  Piezoviscous  flows,  this  assumption  allows  great  simplification  of 
the  Navier-Stokes  equations  and  is  analogous  to  the  result  obtained  by  Denn  [26]  for  cylindrical 
Poiseuille  flow.  However,  we  will  show  that  neglecting  the  cross-film  velocity  component  in  a 
flow  of  a  piezoviscous  fluid  between  parallel  plates  results  in  inconsistent  cross-derivatives  of 
pressure.  To  illustrate,  we  will  begin  by  examining  the  cylindrical  Poiseuille  flow. 
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5.5.1  CYLINDRICAL  POISEUILLE  FLOW 

Denn  [25]  used  a  perturbation  solution  to  establish  that  the  radial  velocity,  v,  was  zero  in 
cylindrical  flow.  In  preparation  for  the  following  section,  we  will  obtain  Denn’s  result  by 
different  means. 

For  isoviscous  liquids,  the  velocity  is  of  the  form 


R2  r2 
u  =  q - q — 

h  2  M  2 


(37) 


and  v  =  0 

where  R  is  the  radius  of  the  cylindrical  duct  and  q  =  4  Q  /  7t  R4 ,  where  Q  is  the  flow  rate. 


The  appropriate  Navier- Stokes  equations  for  an  incompressible  piezoviscous  liquid  are: 

(38) 


dp  du  dp  p  du  d2u 

^TaM¥¥+7a+Me? 


dp  dp  du 

—  =  ap— — 

dr  dx  dr 

where  r  and  £  represent  the  radial  and  axial  directions,  respectively. 
Making  use  of  (37).  equations  (38)  and  (39)  become 

dp  dp 

-  =  -CMqr--2Mq 

dp  dp 

—  =  -apqr— 

CT  Cl 


Combining. 


dp  -2pq 


d%  l-(apqr) 


2  _  2 


dp  2  arp  q 


(39) 


(40) 

(41) 

(42) 

(43) 


dr  l-(a|iqr) 

Taking  the  partial  derivatives  of  (42)  and  (43),  we  obtain  the  following  equations. 

dfdp ^ 


dr 


-  dp  dp  2 
=  2a— 7“  +  ra2 
dE,  dr 


rdpyf(  dp' 


dr 


(37) 


Take  the  partial  derivative  of  (36)  with  respect  to  ^to  obtain: 
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(45) 


cc\cr)  cl  cr  VcriUd 


The  equivalence  of  the  expressions  (44)  and  (45)  for  the  cross  derivative  of  pressure 
ensures  that  the  solution  for  p  is  independent  of  the  integration  path  taken  and  verifies  the 
assumed  velocity  field  given  by  equation  (37).  Therefore  there  is  no  secondary  flow,  v  =  0. 
Note  also  that  in  capillary  flow,  a  pressure  singularity  occurs  where  apqr  =  l  and  that  this 
singularity  will  be  downstream  of  the  location  of  the  unbounded  pressure  associated  with 
piezoviscous  choking  of  flow. 

5.5.2  Flow  Between  Parallel  Plates 

We  will  take  the  same  approach  to  investigate  secondary  flow  between  parallel 
plates.  Here  the  kinematics  can  be  made  more  general  as  the  plates  may  move  relative  to 
one  another  in  the  x-direction.  We  will  assume  a  more  general  form  for  the  velocity  field 
than  the  customary  parabolic  form  of  equation  (28).  Here  we  assume 

u  =  u(z)  and  wsO  (46) 

Remember  v  =  0 ,  since  this  is  a  two-dimensional  problem.  Recall  that  the 
pressure  gradients  are  given  by  (24)  and  (25)  in  general  and  by  (30)  and  (31)  when 
w  =  0 . 


dp  _ 

5x  l-(a|iu2): 


dp  aii1  dp 

dz 


Differentiating  (30)  with  respect  to  z  yields  after  much  manipulation, 

_£_( ££.1  -  3a£E££  |  i  2a2  uz  TgP  Y  & 

dzldxj  dzdx  un  dx  na{dz)  dx 
Similarly,  differentiating  (24)  with  respect  to  x  yields 


(30) 

(31) 


(47) 
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(48) 


-£.(■&).  2a&&  + 


dx  {  dz 


dz  dx 


zx  zzx 


u 


zz  y 


dp  2a  r  dp 
dz  u„ldz. 


dp 

au,—  +  u„ 
2dx  “ 


And  with  ux  =  0,  Equation  (48)  becomes 


—( £El  _  2  a— — +  — T  — 

dxydz)  dz  dx  u  l^dzj  dx 


(49) 


Equating  the  RHS  of  (49)  and  (47)  reveals  properties  of  the  admissible  form  of  u(z)  - 
when  w  =  0 . 

uI22=-aun^  =  -aV'i2uc^  (5°) 

It  is  clear  that  the  customary  parabolic  form  of  the  velocity  is  inappropriate  for  a  *  0, 
since  for  a  polynomial  representation  of  u(z),  a  third  degree  term  is  needed  for  d3u/dz3  to  be 
nonzero. 


Additional  information  can  be  extracted  by  considering  that  du^/dx  =  0,  since  w  =  0. 


Now. 


-a~u2ua 


ap 


^  d2p 


dx1 


=  0 


(51) 


25 


So  that 


a 


(dp)2 
o dx j 


-g2P 

dx2 


Differentiating  eqn.  (30)  with  respect  to  x  gives 

8  (dp 


dx  V5x 

Substituting  (53)  into  (52)  yields 


=  af— T+2a3|i  "2 

IdxJ  % 


3..  u?  I  % 


.dx. 


2  2^P 

Uz2  =  -apuz- 


(52) 


(53) 


(54) 


which  must  be  satisfied  for  whO.  Note  that  the  sign  of  uc  must  be  opposite  that  of  dp/dx . 
This  is  contrary  to  what  is  found  in  practice.  Substituting  the  expression  for  dp/dx  (eqn.  (30)) 
into  (54)  to  obtain 

l-a2jru2  =  -a2p2u2  (55) 


which  is.  of  course,  a  contradiction.  Therefore,  w#0,  and  a  secondary  flow  must  occur  for 
a  *  0  when  a  Poiseuille  component  exists. 


6.  SUMMARY  AND  CONCLUSIONS 

If  care  is  taken  to  avoid  or  correct  for  viscous  heating,  the  results  from  high-pressure 
Couette  rheometers  may  be  used  to  generate  empirical  rate  equations  which  are  useful  in 
modeling  elastohydrodynamic  traction.  This  approach,  however,  fails  to  ellucidate  the  details  of 
lubricant  flow  in  the  concentrated  contact.  The  outstanding  problems  of  surface  roughness 
effects  and  EHD  related  boundary  lubrication  await  analytical  techniques  which  can  handle  the 
change  in  character  of  the  governing  (Navier-Stokes)  equations. 

An  analytical  treatment  of  piezoviscous  liquids  is  presented  which  reveals  that  the 
Reynolds  equation  adequately  captures  the  mechanics  of  a  piezoviscous  fluid  only  when  tcx  « 
1.  In  the  elastohydrodynamic  lubrication  regime,  this  condition  is  satisfied  in  the  inlet  zone 
where  typically  xa  =  0.02.  Consequently,  Reynolds-based  EHD  inlet  analyses  yield  realistic 
predictions  for  film  thickness.  In  the  Hertzian  contact  region,  the  contribution  of  the  pressure 
gradient  across  the  film  may  play  a  very  important  role  and  the  use  of  the  Reynolds  equation  is 
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questionable.  While  the  mathematical  treatment  of  the  subject  becomes  considerably  more 
complex,  these  issues  must  be  dealt  with  if  closure  is  expected  in  predicting  EHD  traction  and 
relating  it  to  independently  determined  lubricant  properties. 
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Figure  2.  Flow  Chart  for  Santotrac  50,  p  «  297  MPa. 
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Figure  3 .  Flow  Chart  for  Mineral  Oil,  LVI 260,  p  =  241  MPa. 


Figure  £  Torque  Transducer  Operating  Principle  for  Rheometer  Shown  in 
Figure  4 
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Figure  6  Flow  Chart  for  a  Polyalphaolefin. 
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Figure  7  a  Flow  Chart  for  H VI 650. 
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Figure  7b .  Flow  Chart  for  H VI 650. 
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Figure  8 .  Flow  Chart  for  Santotrac  40 
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Figure  9  Flow  Chart  for  Grease  and  Base  Oil. 


Figure  10 .  Measured  and  Predicted  Traction  Curves  for  HVI 650. 
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Figure  11.  Mechanical  Shear  Bands  between  Converging  Plates  with  Sliding.  Lower 
Tilted  Surface  is  Stationary. 


variation  of  compressibility  and  pressure  viscosity  coefficient 
pressure. 


